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ABSTRACT 

The concept of (auto) catalytic systems has become a cornerstone in understanding evolutionary processes in var- 
ious fields. The common ground is the observation that for the production of new species/goods/ideas/elements 
etc. the pre-existence of specific other elements is a necessary condition. In previous work some of us showed 
that the dynamics of the catalytic network equation can be understood in terms of topological recurrence rela- 
tions paving a path towards the analytic tract ability of notoriously high dimensional evolution equations. We 
apply this philosophy to studies in socio-physics, bio-diversity and massive events of creation and destruction 
in technological and biological networks. Cascading events, triggered by small exogenous fluctuations, lead to 
dynamics strongly resembling the qualitative picture of Schumpeterian economic evolution. Further we show 
that this new methodology allows to mathematically treat a variant of the threshold voter-model of opinion 
formation on networks. For fixed topology we find distinct phases of mixed opinions and consensus. 
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1. INTRODUCTION 

This paper deals with the understanding of dynamics of systems capable of evolution. In particular we are 
interested in systems of technological evolution with biological evolution as a special case. We are dealing with 
systems of elements (e.g. molecules, species, goods, thoughts) that do not only exist or decay but have the 
possibility to form new elements through recombination of already existing elements. 

For example consider a soup of molecules. Some molecules in the soup can react chemically with each 
other and form new molecules. Some reactions will lead to products that are already existing in the soup. Other 
reactions may produce molecules that are new and therefore transform the soup-environment into something new. 
This in return opens the possibility for chemical reactions that have not been possible before the existence of the 
new product. In this context the concept of the adjacent possible 1 or adjacent probable 2 has been introduced. 
This is the set of objects that can possibly become produced within a given time span into the future. What 
can be produced in the next timestep depends crucially on the details of the set of elements that exist at a given 
time. Existing elements undergo an evaluation within the newly formed context. 2 

In this paper we are basically interested in three types of phenomena. First, cascades of creation, second, 
cascades of destruction. The existence of these kind of dynamics has been postulated in the Schumpeterian 
view on economic evolution 3 long time ago. Third, we study formation processes in evolutionary dynamics. In 
this scenario elements exist in two different states 0/1 (like existing/not existing, on/off, active/inactive). Each 
element is connected to a set of other elements through a (fixed) interaction matrix. An element can change its 
state depending on the states of the elements it is connected to. For example, some type of molecule may change 
its state from non existing to existing when sufficiently many other molecules also exist - or vice versa. This 
problem can be cast onto an opinion formation 4 10 scenario. 11 

One of the substantial problems in dealing with evolutionary systems is the notoriously high or even unbound 
dimensionality of systems undergoing evolution. Just think of how many biological species exist on this planet 
or how many different goods and services are available, e.g. in Europe. Every element (good/species/chemical 
substance) i is present in an evolving system with some relative frequency X{. X{ is the normalized abundance of 
the element z, such that x % ~ 1- Think of each X{ as an entry in a d- dimensional vector containing all possible 



elements in the universe. Inter-dependence between pairs of elements can be coded in a large interaction matrix 
a. If, for instance, i uses a service provided by j we set — 1 otherwise = 0. For ternary (and higher) 
relations a is not a matrix but a tensor which we call the production rule table or simply rule table. Suppose we 
look at chemical reactions then = 1 when i is the product of a chemical reaction of the substrates j and 
and ctijk = if this chemical process is not permitted by the rules of chemistry. 

High dimensionality usually makes such systems particularly difficult to treat in an analytic fashion. On the 
other hand the combinatorial explosion of elements that emerges through the recombinations of other elements, 
severely limits computational analysis as well. The exact way one element influences another is often only poorly 
known, or not known at all. For example, we do not exactly know how many biological species are living on this 
planet right now (although we may make educated guesses). Needless to say that we know even less about the 
actual interactions and relations that species involve in. Consequently we usually can not exactly describe the 
inter-dependencies between species (e.g. trophic dependencies) which govern their dynamics. We are certainly 
incapable of prestating future inter-dependencies of species that might come into being in the future course of 
earth history. This is to say, that we usually know very little about the true structure of the interaction matrices 
or tensors a. However, we may happen to know several things, constraining our ignorance. For example, we 
might have knowledge about the average number of chemical reactions producing the same chemical substance, 
or the average number of trophic dependencies an individual in an ecosystem is maintaining. 

It was suggested 12 to meet the above problems of dimensionality and ignorance by two model assumptions: 
(i) The focus is shifted from the actual frequency of elements to the system's diversity, where diversity is 
defined as the number of existing elements. An element exists, if X{ > 0, and does not exist if X{ = 0. (ii) The 
entries of the interaction matrices or rule tables a that encode the relation between the specific elements of the 
system, are modelled as random. For simplicity the entries of a are or 1. 

These two assumptions are powerful tools in mapping the dynamics - that ideally should be described on a 
level of differential equations of the relative frequencies (concentrations) of elements - onto recurrence relations 
for the system diversity or other variables of interest. 

In the following we present two examples of how evolution dynamics can be treated methodologically in this 
manner. The first deals with cascades of creation and destruction in evolutionary systems. This is done in 
section [2] where we also show how to extend the analysis of earlier work. 12,13 The second example deals with the 
formation process in section [3J 



We have analyzed systems where pairs of elements exist that can be combined to form a third element. 12, 13 
Such systems are sometimes called catalytic networks. We have demonstrated the existence of both creative 
and destructive cascades (algorithmic phase transitions). We first review the results and then extend them 
to the situation where n elements can be combined to form a product-element. Note, that this model is a 
mathematically tractable variant of the so called bit-string model of biological evolution. 1 

Suppose there is a maximal diversity of elements in the system d that can possibly exist and write the relative 
frequencies xi (how much there is of i) of the elements i as the d dimensional vector x = (xi,X2, • • • , This 
vector x describes the state of the system. The dynamics of x is governed by the famous catalytic network 
equation 



where the second term ensures normalization of x. The tensor elements serve as the 'rule table', telling 
which combination of two elements j and k can produce a third (new, different) element i. The rule table element 
dijk is the rate at which element i can get produced, given the elements j and k are abundant at their respective 
concentrations Xj and x^. We call a pair of elements (j, k) productive when there exists an i such that > 
is strictly larger than zero. Special cases of Eq. ([I]) are the Lotka Volterra replicators, 14 the hypercycle, 15 or the 
Turing gas. 16 Equation (p]) has been analyzed numerically, 17, 18 however for rather moderate system sizes. 



2. CREATIVE AND DESTRUCTIVE CASCADES 




(1) 



Many goods and services require more than two inputs. For instance, a car-engine consists of many parts that 
have to be put together to form a functional entity, or the production of wine requires barrels, grapes, bottles, 
human labor and appropriate quantities of sun and rain. Therefore we generalize the catalytic network equation 
to processes, where n elements instead of two elements are necessary to form a new one. To do so, we have to 
pass from considering index pairs (j, k) like in Eq. ([!]), to n dimensional index vectors j = (ji, . . . , j n ) to indicate 
the necessary n substrate elements. Therefore the rule table becomes a tensor = &i,j 1 ,j 2 ,---,jn depending on 
the product i and n substrate-elements j m , where m = 1,2, . . . , n. To model the effect of the availability of 
elements on the production we also have to pass from the frequency pair XjXk in Eq. ([!]), to the product of the 
n frequencies Xj ± Xj 2 . . . Xj n . The notion of the productive pair naturally extends to the notion of productive sets 
such that a set of elements j is productive when there exists an i such that > 0. Replacing this appropriately 
into Eq. gives us the generalized catalytic network equation. Remarkably, this generalization does not alter 
the arguments that allow us to analyze the creative and destructive cascades of production we are interested in. 

2.1 Cascades of creation 

First of all, applying (i) from the previous section, we reduce the dimensionality of the problem by counting the 
elements X{ > 0, and define the relative diversity a t G [0, 1] of the system such that a t d is the number of existing 
elements at time t. We assume that the system has a growing mode only (tensor elements a > 0, i.e. a = 1 
but never negative). Entries 1 in the rule-table tensor a indicate the productive sets. Suppose, applying (ii), we 
only know the number r such that rd is the number of such productive sets occurring in the system and that 
the entries 1 and in a are distributed completely randomly, r constitutes our constraining knowledge on a. 
The only variables characterizing our system therefore are r and n. Let us suppose - as the initial condition - 
that initially, at time t = 0, only a few elements ao d exist. This means that a^d elements in the initial element 
vector are different from zero. We are interested in the asymptotic value (t — > oo) of the relative diversity a^. 

The fraction of sets of n elements in the system which are productive is given by r d divided by the binomial 

coefficient of d over n elements ^ ^ j . In a given set of size a t d there are ^ a ^ ^ sets consisting of n elements 

to be found on average. Due to the assumed randomness of the distribution of productive sets of elements the 
number of produced elements is 

Yet, about r d a™_ 1 of those elements already have been produced in the preceding time step t — 1 and are certainly 
contained in a t already. They consequently have to be subtracted as already existing. Due to the randomness 
of the production of elements only a fraction 1 — a t of the remaining elements will not be present in a t and truly 
be newly created. Putting all this together yields 

a £+ i=a £ + Aa £ , Aa t = r (1 - a t ) (a? - a?_ x ) . (3) 

at timestep t = we have to set a_i =0, by convention since ao is the initial condition. 

These equations are exactly solvable in the long-time limit. For this end define, c t = Aa^+i/Aa^, and look 
at the asymptotic behavior, c = lim^oo c t . From Eq. Q we get 

c = nr(l-a 00 )a 1 ^ 1 . (4) 

On the other hand we can estimate asymptotically by 

CO 

Ooo = a c l = —5- . (5) 

1 — C 

£=0 

By inserting Eq. Q into Eq. ([5| one gets a N + l'th order equation 

«oc - (M) = nr(l - a 00 )a^ , (6) 



Figure 1. Phase diagram of the creative dynamics over the r-ao plane for n = 2 and n = 3. The diagram shows clearly 
how the system starts to exhibit a phase transition for r ~ 2 in the n — 2 case and how the critical initial complexity ao 
where the phase transition occurs gets smaller for increasing values of r. The phase transitions occurs for larger values of 
ao in the n = 3 case. 



whose solutions are the solution to the problem. Most remarkably these solutions are mathematically identical 
to the description of real gases, i.e. Van der Waals gases, independent of the parameter n. One can make the 
relation to the Van der Waals gas explicit by defining, V = a^, T = ao(nra^+ 1 ) -1 , and P = a^ 1 + (nra^ 1 ) -1 . 
This allows to generically map Eq. (|6| onto the famous equation 

( p -^) V = T • (7) 

which is exactly a Van der Waals gas of point-particles with constant internal pressure. The meaning of 'pressure' 
and 'temperature' in our context is intuitively clear. As comparable to real gases our n-ary assembly system 
shows a phase transition phenomenon which are of the same universality class independent of the choice of 
n. The corresponding phase diagram, as a function of the model parameter r and the initial condition ao is 
shown in Fig. [I] for the choice n = 2 and 3. Note that assembly systems with larger values of n undergo the 
phase transition at larger values of ao for fixed r so that larger sets of random entities are required to trigger 
self-sustaining creative cascades. 

The update-equations (|3| are scale invariant. The solutions of the problem in terms of relative diversities do 
not explicitly depend on the upper bound of diversity d and the results therefore apply also in the limit d — > oo. 
One can think of open systems which can not be described as a limit of bounded systems. Such systems can be 
treated in a similar fashion, however this task is of considerably higher mathematical complexity. Work dealing 
with this problem is in preparation. 

2.2 Cascades of destruction 

In the dynamics studied so far, diversity can only increase due to the positivity of the elements in a. It is 
important to note that in this setting the phase transition can not be crossed in the backward direction. This 
is because of two reasons. First, the system forgets its initial condition ao once it has reached the (almost) fully 
populated state. This means that after everything has been produced one can not lower the initial set size any 
more. In terms of the Van der Waals gas equation analogy we can not lower the 'temperature'. Second, if r and 
n are homogeneous characteristics of the system then it is also impossible to manipulate the 'pressure' of the 
system and we remain in the fully populated phase for ever. The natural question thus arises what happens to 
the dynamics if one randomly kills a fraction of elements in the fully (or almost fully) populated phase. In the 
case that an element k gets produced by a single productive set j and one of the members of the j m gets killed, k 
can not be produced any longer. We call the random removal of i a primary defect, the result - here the stopped 
production of k - is a secondary defect. The question is whether there exist critical values for r and the primary 
defect density Oq, such that cascading defects will occur. 



As before one may approach this question iteratively, by asking how many secondary defects will be caused 
by an initial set of Sod randomly removed elements in the fully populated phase. The possibility for a secondary 
defect happening to element k requires that all productive sets, which can produce fc, have to be destroyed, i.e. 
at least one element j m of the productive set has to be eliminated for each productive set j. On average there 
are r such productive sets for each element i. Let us introduce the first defect. This defect will on average affect 
nr productive sets in the system, i.e., there will be nr elements that lose one way of being produced. This is 
so, since there are dr productive sets and there are ndr indices referring to an element involved in denoting 
the productive set. Consequently there are nr indices on average per element. To look how adding primary 
defects effects the system and the number of elements that have lost one, two, and more ways of getting produced 
iteratively requires some 'book- keeping'. However, from here on we basically merely have to adapt a derivation 
presented in earlier work 13 in a straight forward fashion. Since the mathematics is not difficult but lengthy we 
directly jump to the result 

Si=7(r)f rn (S )5 r , (8) 

with 

7(0 = ^(r-ljl-V-- 1 

(9) 

frn(So) = ^O^t-^of • 

The result remains to be iterated from secondary to tertiary and higher order defects. To do so we think of 
collecting the primary and secondary defects together and assume that we would start with a new primary defect 
set of relative size S' = So + Si . The tertiary defects therefore gets estimated by S2 = S[ — Si , where S[ are the 
secondary defects associated with S' . This leads to the recursive scheme, 

m 

A m = J2 s k , S m+ i = 7 (r)/ rn (A m )A^ - A m + (5 , (10) 

k=0 

where A t are the culminated defects after t timesteps. The result in terms of a phase diagram of the iteration 
scheme is given in Fig. [2] for n = 2 and 3. The asymptotic culminated total defect size Aoo (for t — > 00) is 
shown as a function of the parameters r and the initial defect density So. As before a clear phase transition is 
visible, meaning that at a fixed value of r there exists a critical number of initial defects at which the system 
will experience a catastrophic decline of diversity. 

It is interesting that for the complete destruction of diversity (plateau in Fig. |2| not very large values of So 
are necessary. We also want to stress that assembly systems with larger values of n collapse at smaller values of 
So- They are more sensitive to the exposure to initial defects. 

In a more realistic scenario one can assume that (exogenous) random processes are introducing new elements 
at a given rate into a system just as other (exogenous) random processes may introduce defects. We have shown 13 
how, being equipped with the update equations for the creative and the destructive cases Eqs. ([3| and (|8|, these 
update equations can be coupled to study the combined dynamics driven by exogenous random events. This 
allows to take a look at how diversity of systems may evolve over time, driven by the spontaneous creation and 
destruction processes, which may reflect exogenous influences. There we have also shown that the increment 
distribution of Aa t = a t — a t -i is power-law distributed, when the system is prepared critically, i.e. at the onset 
of the creative phase-transition. 

We may conclude from our analysis of generalized catalytic systems that average number of elements required 
to build a new element of the system does not change the overall picture previously drawn 12, 13 although larger 
n make it more difficult to reach the creative phase transition in the low diversity phase and easier to reach the 
destructive phase-transition in the high diversity phase. The limit n — » 1 is comparable with the situation of 
the opinion formation problem under the unanimity-rule 11 for an average connectivity r of each node. Since the 
stability of biological networks are often analyzed under the unanimity-rule governing trophic dependencies, 19 21 
where a species can survive as long as it has one remaining trophic link, it is interesting to see that this is 
the most stable model-assumption one can pose (n = 1); this is to say trophic-dependencies-only is a friendly 
scenario. If the survival of species are also depending on other than trophic factors this can only result in an 
increase of the parameter n to values larger than one. Thus biological networks can be considerably less stable 
than the analysis, based on trophic dependencies only, would suggest. 



n=2 



n=3 
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Figure 2. Phase diagram of the destructive dynamics over the r-5o plane for n = 2 and n — 3. We see that relative 
small numbers of initial defects suffice for the whole system to break down. The culminated defects Aoo first raise linearly 
dominated by the number of primary defects 5o . At the critical of primary defects the system breaks down almost without 
precursor in the number of culminated defects. We may also see how the n = 3 case is considerably more sensible to 
defects than the n — 2 case. 



3. FORMATION DYNAMICS 

So far we have been interested in creation and destruction processes on catalytic sets. We now turn to the analysis 
of the formation processes in evolutionary dynamics. It was shown that the formation problem of systems of 
evolution can be mapped onto an opinion- formation problem. 11 As mentioned in the introduction, elements in 
this scenario exist in two different states 0/1 (e.g. existing/not existing) and each element is connected with a set 
of other elements through a (fixed) interaction matrix. An element can change its state depending on the states 
of the elements it is connected to. For example, some type of molecule may change its state from non existing to 
existing when sufficiently many other molecules, from which it can be produced in chemical reactions, also exist - 
or vice versa. Another example, which can be understood as an opinion-formation problem, is the ongoing debate 
about the standard for high definition DVDs (Blu-ray against HD-DVD), paralleling the historical Betamax vs. 
VHS conflict. We assume that one is more likely to favor one alternative if this alternative is also favored in one's 
local neighborhood. If all of my friends would have shared their tapes in the VHS standard (because a certain 
branch of movies is only released in this format), it is more likely that I would have dumped my Betamax video 
recorder and used VHS from then on. It is also possible that two different opinions can co-exist, e.g. Mac and 
Windows users. In the next section we will propose a model capable of showing these phenomena. A node will 
change its internal state if the fraction of neighbors holding the opposite opinion (state) is above a predefined 
threshold 0.5 < h < 1. These opinion formation models are sometimes called threshold voter models. 22,23 




Figure 3. A pictorial description of the update process with h — 0.8. The node marked with an arrow gets updated, it 
has five neighbors for which we show two possible configurations of internal states, (a) Three out of five neighbors are in 
a different state, the node stays unchanged, (b) Four out of five neighbors are in a different state leading to a change in 
the internal state of the updated node. 
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Figure 4. The order parameter cij ry-j as a function of the initial population sizes, a[], for N = 10 4 , ft = 0.8. (a) k = 2 
for all nodes (unanimity regime), (b) random network with k = 10 (intermediate regime) and (c) random network with 
k = 9x 10 3 . With increasing k the transition between the phase of coexistence and consensus becomes sharper. 



3.1 The model 

We represent an agent i as a node in a network composed of N agents. The internal state of the node corresponds 
to the agent's choice, e.g. which technology he uses in every-day situations. For simplicity we break down each 
decision into binary choices, yes/no, VHS/Betamax, MS Windows /Mac, etc. Connected nodes influence each 
other in the decision-making process, i.e. they tend to persuade their neighbors to adopt their own choice. The 
decision-making process consists of four steps (see Figj3|: 

• Pick a node i randomly. 

• Check the internal state of each of z's neighbors. 

• If the fraction of states holding opinion 1(0) is above the threshold ft, i adopts opinion 1(0). 

• If the update threshold is not exceeded, i keeps its initial state. 



We assume that the only information available about the topology of the network is the probability for any 
two nodes to be connected. Our ignorance thus tells us, following model assumption (ii), to model the interaction 
matrix as a random network 24 where N nodes are randomly linked with L links (self-interactions are forbidden). 
Of course any other network topology could be used, however at the cost of more mathematical trouble. 25 If two 
nodes i and j are connected, we set the corresponding entries (i, j) and (j, i) to 1, if two nodes do not interact 
the entry would be set to 0. For sufficiently large systems the crucial parameter is the average connectivity 
k = L/N and our results are independent of the actual number of agents. The update threshold necessary for a 
node's change of opinion, ft, has to be higher than 0.5 in order to be meaningful in the above sense. The update 
is carried out random sequentially. In a network containing N agents at any time t, there are A® nodes with 
choice and A] nodes with choice 1. We again use dimensionless variables and define the fraction of nodes in 
state 0/1 as a^ 1 = A^ 1 /N. The role of the a^ 1 is comparable to the role of the relative diversity a t in the 
previous section, and reduce the complexity of the problem according to model assumption (i). If the sequence 
of updates is complete we increase t by 1 and start the next random sequence of updates. We are interested in 
the order parameter lim^oo a^ 1 — aj^. 

3.2 Theory and simulations 

We now derive a master equation for the evolution of this system by adopting the same philosophy as before. 
First, let us outline the main difference with respect to the previous scenario. In the creative dynamics we 
incorporated updates in one direction, e.g. from not-produced — > produced. In the formation dynamics the 
updates are required to take place in both directions. Thus our aim is to derive the transition probabilities 
A^ 1 and Aj - * for both processes, an update from to 1 and vice versa with combinatorial arguments. 12 We 
start at t = where we have a fraction of a® nodes in state 0. The probability that at time t an update from 
to 1 occurs is denoted by Aq^ 1 . This transition rate is the sum over the probabilities for all configurations of 



neighbors where at least a fraction of h neighbors is in state 1 and we simultaneously have picked a node initially 
in state 0, 

K^=a°o £ (fWl-agyK)^ , (11) 



i=\kh\ 



where |\] denotes the ceiling function, the nearest integer being greater or equal. To derive Aq one just has 



to exchange and 1 in Eq.(ll ). The equation governing the system's evolution from t = to t = 1 thus is given 

4 = a° + Aj-° - A"- 1 . (12) 

Let us examine some special cases now. 
3.2.1 Unanimity regime 

If k and the update threshold h are chosen such that \kh] = k holds, the evolution halts after one iteration. 
We call this the unanimity regime. Here the update rule boils down to the unanimity rule, i.e. all linked nodes 
have to be in the same internal state to allow an update. Let the following account for an argument: Consider 
the second iteration. Assume that no consensus has been reached and we thus can find two neighboring nodes 
favoring different choices, say agent i is in state 0, j holds state 1. i could be updated if each of his other 
neighbors is in state 1. But if this was the case, i could not be in 0. Either the update to 1 would have occurred 
in the first iteration, or there is an agent k in z's neighborhood which also holds state and can not be updated 
due to his connectedness with i. So, as one can easily check by inspection, the dynamics of the system freezes 



after the first iteration. For the special case k = 2 (constant) we find that = a\. Inserting this in Eq.(12) 
yields 

«°oo = 3(a°) 2 - 2(a°) 3 . (13) 



Results of a simulation on a regular ID circle network with N = 10 4 confirm the theoretical prediction of Eq.( 13), 
see FigJ^a). 

3.2.2 High connectivities 

For higher values of k and h the gridlock scenario described above does not appear. In subsequent iterations we 
still observe updates. However, after waiting sufficiently long the evolution converges to a halt. As the fraction 
of nodes in one state increases, this raises the probability that the respective fraction of nodes will grow even 
further (as the words spread amongst the agents). This may force the system into consensus. Let us consider 
arbitrary time steps now. To derive A^ - * 1 we look at the number of updates being possible with respect to the 
starting conditions, and we subtract the probabilities that these have already occurred at former time steps. 
This directly leads to 

A?- 1 = a° ( ± (J) (1 - 4Y (a?)*"' - t (J) (1 - aUY (aUfA ■ (14) 
\i=ffe/i] V J i=\kh~\ V J j 

Without being able to offer any more surprises we can finally state the master equation for arbitrary times, 
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In Fig|4|c) we test Eq.(15) against simulations for a system with N = 10 4 , h = 0.8 and k — 9 x 10 3 and find a 
formidable level of agreement. The plot allows us to identify three regimes. There are two regions of consensus 
with a sharply separated linear regime of coexistence. 

For the case of a fully connected network we can state the asymptotic population sizes explicitly. If the initial 
populations are both smaller than the update threshold, i.e. <2q < h and 1 — a[j < /i, there will be no updates at 
all. Otherwise the system reaches consensus. So, for k = N — 1 we find 

a °oc = a °o + (1 - a °o) © Oo - h ) ~ «o@ (1 - a °o ~ h ) , (16) 



(a) (b) 




Figure 5. as a function of the initial fraction size a[] and connectedness k/N . Simulations were performed on random 
networks with N = 10 3 and h = 0.8. Two symmetrical regions of consensus and a mixed phase in between are observed, 
(a) Phase diagram for the unanimity and the intermediate regime. With increasing connectedness the regions become 
more distinguishable, (b) For even higher values of the connectedness the dependence of the initial fraction size becomes 
linear in the mixed phase. 



where 9 (x) is the Heaviside step function. How can we understand the transition between this curve and the 
smooth solution found in FigQa)? 

Intermediate regime. An interesting case appears when we leave the unanimity regime but keep the connec- 
tivity low compared to N. We have performed simulations with N = 10 4 , h = 0.8 and k — 10. Results are shown 
in FigQb) and compared to the theoretical prediction of Eq.(15). Although the agreement can be coined fair 
on a sunny day, it is not as good as in the previous cases. The question arises what is missing in our description. 
In the iterative scheme we effectively calculate the relative numbers of updates without assuming any kind of 
correlations between the configurations at two different time steps. When we know that an update has taken 
place we have revealed some information about the neighborhood of the updated node. As an extremal case 
of this we have seen in the unanimity regime that after the first iteration no more updates can take place. In 
our analytical work we have not included this knowledge and assume the local neighborhoods to be maximally 
random with respect to the constraining population sizes. While this is a fair assumption for the vast majority 
of values for fc, there exists an intermediate regime where this is not the case. The dynamics of the system is 
shown in the phase diagram, Figj5[ It illustrates the size of the respective regimes and their dependence on the 
parameters a® and connectedness k/N. 

For the voter-model it is known that the system will always end up in a consensus state. In contrast to this 
we also find stable solutions for coexistence. This sheds some light on the question of which parameters are 
controlling whether a collective decision-making process ends in consens (like in the VHS-Betamax example) or 
in a coexistence (e.g. the Mac and Windows users). 

This is possible due to two ingredients in the model presented here. First, our updates are carried out random 
sequentially. As a consequence the distribution of opinions among the nodes does not stay maximally random 
(w.r.t. initial conditions) and groups of agents sharing the same opinion appear. Second, the threshold h can 
also be viewed as a level of tolerance in the sense of allowing a node to maintain its internal state even though 
a majority of the neighbors thinks different. 



4. CONCLUSIONS 

Common ground in systems capable of evolution is that for the production of new species/goods/ideas etc. the 
existence of specific other items is a necessary precondition. In this context (auto) catalytic systems has become 
a cornerstone in understanding evolutionary processes in various fields. Notoriously high dimensionality of 
problems and poor knowledge on the exact form of interactions governing such systems complicate the situation. 
We have shown in previous work how reducing the description of dynamics from the detailed level of the frequency 
of elements to the level of system diversity and how modeling interactions as random (when no better knowledge 



is available) paves a path towards the analytic treatability of systems of evolution. We present two examples 
to make clear the power of the arguments and the wide applicability of the method. First, we applied the 
philosophy to massive events of creation and destruction in technological and biological networks. Here we 
extended and generalized earlier results on catalytic systems. We find cascading events. Generalized catalytic 
systems can pass from low but critically large initial diversity to high diversity and they can break down from 
high to low diversity at a critical point of primary defects (removal of elements). The existence of this kind of 
phenomena in evolution is long known 26 in biology (e.g. the Cambrian explosion) and has been formulated in 
heuristic and qualitative terms in the economical context by J. A. Schumpeter. With our work we increase our 
ability to describe his 'creative gales of destruction' in a more quantitative framework. Further we show that the 
new methodology allows to mathematically treat a variant of the threshold voter-model of opinion formation on 
random graphs. For fixed topology, we find distinct phases of mixed opinions and consensus. Stated in terms of 
systems of evolution the phases of consensus correspond to the phases of maximal and minimal diversity where 
either (almost) nothing or (almost) everything possible exists. The stable mixed opinions correspond to a phase 
where only a fraction of what possibly can be produced, exists in a stable way. Whereas for the voter dynamics 
it is known that it always leads to a consensus state (i.e. only everything/nothing or one of two opinions/states), 
the model studied here allows the coexistence of two states. 
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